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Abstract 

The paper develops DILOC, a distributive, iterative algorithm that locates M sensors in W^,m > 1, with 
respect to a minimal number of m + 1 anchors with known locations. The sensors exchange data with their neighbors 
only; no centralized data processing or communication occurs, nor is there centralized knowledge about the sensors' 
locations. DILOC uses the barycentric coordinates of a sensor with respect to its neighbors that are computed 
using the Cayley-Menger determinants. These are the determinants of matrices of inter-sensor distances. We show 
convergence of DILOC by associating with it an absorbing Markov chain whose absorbing states are the anchors. 
We introduce a stochastic approximation version extending DILOC to random environments when the knowledge 
about the intercommunications among sensors and the inter-sensor distances are noisy, and the communication links 
among neighbors fail at random times. We show a.s. convergence of the modified DILOC and characterize the error 
between the final estimates and the true values of the sensors' locations. Numerical studies illustrate DILOC under 
a variety of deterministic and random operating conditions. 
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I. Introduction 

Localization is a fundamental problem in sensor networks. Information about the location of the sensors is key to 
process the sensors' measurements accurately. In applications where sensors are deployed randomly, they have no 
knowledge of their exact locations, but equipping each of them with a locaUzation device like a GPS is expensive, 
not robust to jamming in military appUcations, and is usually of hmited use in indoor environments. Our goal is to 
develop a distributed (decentralized) localization algorithm where the sensors find their locations under a hmited 
set of assumptions and conditions. To be more specific, we are motivated by apphcations where N = M + m + l 
sensors in (for example, ro = 2 corresponds to sensors lying on a plane, while for m = 3 the sensors are in 
three dimensional Euclidean space) are deployed in a large geographical region. We assume that the deployment 
region hes in the convex hull of a small, in fact minimal, number of m + 1 anchors, (m + 1) <C M. The anchors 
know their locations. In such situations, the large geographical distances to the anchors makes it highly impractical 
for the M non-anchor sensors to communicate directly with the anchors. Further, to compute the locations of the 
non-anchor sensors at a central station is not feasible when M is large, as it requires large communication effort, 
expensive large-scale computation, and adds latency and bottienecks to the network operation. These networks call 
for efficient distributed algorithms where each sensor communicates directiy only with a few neighboring nodes 
(either sensors or anchors) and a low order computation is performed locally at the sensor and at each iteration 
of the algorithm, for example, see [1]. In this paper, we present a Distributed Iterative Localization algorithm 
(DILOC, pronounced die-lock) that overcomes the above challenges in large-scale randomly deployed networks. 

In DILOC, the sensors start with an initial estimate of their locations, for example, a random guess. This random 
guess is arbitrary and does not need to place the sensors in the convex hull of the anchors. The sensors then 
update their locations, which we call the state of the network, by exchanging their state information only with a 
carefully chosen subset of m -|- 1 of their neighbors. This state updating is a convex combination of the states of 
the neighboring nodes. The coefficients of the convex combination are the barycentric coordinates [2], [3] and are 
determined from the mutual inter-sensor distances among the sensors using the Cayley-Menger determinants. At 
each sensor I its neighborhood set contains m+1 sensors, for example, its closest m+1 sensors, such that sensor I 
lies in the convex hull of these m+1 neighbors. These neighbors may or may not include the anchors. 

DILOC is distributed and iterative; each sensor updates locally its own state and then sends its state information 
to its neighbors; nowhere does DILOC need a fusion center or global communication. We prove almost sure (a.s.) 
convergence of DILOC in both deterministic and random network environments by showing that DILOC behaves 
as an absorbing Markov chain, where the anchors are the absorbing states. We prove convergence under a broad 
characterization of noise. In particular, we consider three types of randomness, acting simultaneously. These model 
many practical random sensing and conmiunication distortions as, for example, when: (i) the inter-sensor distances 
are known up to random errors, which is common in cluttered environments and also in ad-hoc environments, where 
cheap low resolution sensors are deployed; (ii) the communication hnks between the sensors fail at random times. 
This is mainly motivated by wireless digital conmiunication, where packets may get dropped randomly at each 
iteration, particularly, if the sensors are power limited or there are bandwidth or rate conmiunication constraints in 
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the network; and (iii) the communication among two sensors, when their communication link is active, is corrupted 
by noise. 

Although a sensor can only communicate directly with its neighbors (e.g., sensors within a small radius), we 
assume that, when the Unks are deterministic and never fail, the network graph is connected, i.e., there is a 
communication path (by multihop) between any arbitrary pair of sensors. In a random environment, inter-sensor 
communication links may not stay active all the time and are subject to random failures. Consequently, there may 
be iterations when the network is not connected; actually, there might never be iterations when the network is 
connected. We will show under broad conditions almost sure convergence of an extended version of DILOC that 
we term as Distributed LocaUzation in Random Environments (DLRE). DLRE employs stochastic approximation 
techniques using a decreasing weight sequence in the iterations. 

In the following, we contrast our work with the existing Uterature on sensor localization. 

Brief review of the literature: The literature on localization algorithms may be broadly characterized into 
centraUzed and distributed algorithms. Illustrative centraUzed locaUzation algorithms include: maximum Ukelihood 
estimators that are formulated when the data is known to be described by a statistical model, [4], [5]; multi- 
dimensional scaUng (MDS) algorithms that formulate the localization problem as a least squares problem at a 
centraUzed location, [6], [7]; work that exploits the geometry of the EucUdean space, Uke when locating a single 
robot using trilateration in m = 3— dimensional space, see [8] where a geometric interpretation is given to the 
traditional algebraic distance constraint equations; locaUzation algorithms with imprecise distance information, 
see [9] where the authors exploit the geometric relations among the distances in the optimization procedure; for 
additional work, see, e.g., [10], [11]. Centralized algorithms are fine in small or tethered network environments; 
but in large untethered networks, they incur high communication cost and may not be scalable; they depend on the 
availability and robustness of a central processor and have a single point of failure. 

Distributed localization algorithms can be characterized into two classes: multilateration and successive refine- 
ments. In multilateration algorithms, [12], [13], [14], [15], each sensor estimates its range from the anchors and 
then calculates its location via multilateration, [16]. The multilateration scheme requires a high density of anchors, 
which is a practical limitation in large sensor networks. Further, the location estimates obtained from multilateration 
schemes are subject to large errors because the estimated sensor-anchor distance in large networks, where the anchors 
are far apart, is noisy. To overcome this problem, a high density of anchors is required. We, on the other hand, 
do not estimate distances to far-away nodes. Only local distances to nearby nodes are estimated that have better 
accuracy. This allows us to employ a minimal number of anchors. 

A distributed multidimensional scaling algorithm is presented in [17]. Successive refinement algorithms that 
perform an iterative minimization of a cost function are presented in [18], [19], [20]. Reference [18] discusses 
an iterative scheme where they assume 5% of the nodes as anchors. Reference [20] discusses a Self-Positioning 
Algorithm (SR\) that provides a GPS-free positioning and builds a relative coordinate system. 

Another formulation to solve localization problems in a distributed fashion is the probabiUstic approach. Non- 
parametric belief propagation on graphical models is used in [21]. Sequential Monte Carlo methods for mobile 
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localization are considered in [22]. Particle filtering methods have been addressed in [23] where each sensor stores 
representative particles for its location that are weighted according to their likelihood. Reference [24] tracks and 
locates mobile robots using such probabilistic methods. 

Completion of partially specified distance matrices is considered in [25], [26]. The approach is relevant when the 
(entire) partially specified distance matrix is available at a central location. The algorithms complete the unspecified 
distances under the geometrical constraints of the underlying network. The key point to note in our work is that 
our algorithm is distributed. In particular, it does not require a centralized location to perform the computations. 

In comparison with these references, our algorithm, DILOC, is equivalent to solving in a distributed and iterative 
fashion a large system of linear algebraic equations where the system matrix is highly sparse. Our method exploits 
the structure of this matrix, which results from the topology of the coimnunication graph of the network. We prove 
the a.s. convergence of the algorithm under broad noise conditions and characterize the bias and mean square error 
properties of the estimates of the sensor locations obtained by DILOC. 

We divide the rest of the paper into two parts. The first part of the paper is concerned with the deterministic 
formulation of the locahzation problem and consists of sections II-IV. Section II presents preUminaries and 
then DILOC, the distributed iterative localization algorithm, that is based on barycentric coordinates, generaUzed 
volumes, and Cayley-Menger determinants. Section III proves DILOC's convergence. Section IV presents the 
DILOC-REL, DILOC with relaxation, and proves that it asymptotically reduces to the deterministic case without 
relaxation. The second part of the paper consists of sections V-VI and considers distributed localization in random 
noisy environments. Section V characterizes the random noisy envirormients and the iterative algorithm for these 
conditions. Section VI proves the convergence of the distributed locahzation algorithm in the noisy case that relies on 
a result on the convergence of Markov processes. Finally, we present detailed numerical simulations in Section VII 
and conclude the paper in Section Vni. Appendices I-III provide a necessary test, the Cayley-Menger determinant 
and background material on absorbing Markov chains. 

II. Distributed Sensor Localization: DILOC 

In this section, we formally state DILOC (distributed iterative locahzation algorithm) in m-dimension Euclidean 
space, (m > 1), and introduce the necessary notation. Of course, for sensor localization, m = 1 (sensors in a 
straight fine), m = 2 (plane), or m = 3 (3d-space.) The generic case of m > 3 is of interest, for example, when 
the graph nodes represent m-dimensional feature vectors in classification problems, and the goal is still to find in 
a distributed fashion their global coordinates (with respect to a reference frame.) Since our results are general, we 
deal with m-dimensional 'locahzation,' but, for easier reference, the reader may consider m = 2 or m = 3. To 
provide a quantitative assessment on some of the assumptions underlying DILOC, we will, when needed, assume 
that the deployment of the sensors in a given region follows a Poisson distribution. This random deployment is 
often assumed and is reahstic; we use it to derive probabihstic bounds on the deployment density of the sensors 
and on the conmiunication radius at each sensor; these can be straight forwardly related to the values of network 
field parameters (hke transmitting power or signal-to-noise ratio) in order to implement DILOC. We discuss the 
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computation/communication complexity of the algorithm and provide a simphstic, yet insightful, example that 
illustrates DILOC. 

A. Preliminaries and Notation 

Recall that the sensors are in M™. Let 6 be the set of sensors or nodes in the network decomposed as 

e = KU^], (1) 

where k is the set of anchors, i.e., the sensors whose locations are known, and f2 is the set of sensors whose 
locations are to be determined. By | • | we mean the cardinality of the set, and we let |6| = N, = m + 1, 
and |fi| = M. For a set ^ of sensors, we denote its convex huU by C {^Y- For example, if is a set of three 
non-coplanar sensors in a plane, then C {^) is a triangle. We now define a few additional sets needed. 

Let dik be the Euclidean distance between two sensors l,k G Q. We associate with the sensor I G fl, a positive 
real number, r; > 0, and two sets /C {I, ri) and 6; (ri): 

lC{l,n) = {keQ: dik<ri}, (2) 

e;(n) c/c(Z,n), l^Oiin), /eC(e;(n)), 16; = m + 1, ^e,(n)^o> (3) 

where ^e,(r,) is the generalized volume (area in rn = 2-d, volume in to = 3, and their generalization in higher 
dimensions) of C (G; (n)). The set IC{l,ri) groups the neighboring sensors of / within a radius ri and, by (3), 
0; (fi), which we will often represent simply as O;, assuming r; is understood from the context, is a subset of m,+ 1 
sensors such that sensor I lies in its convex hull but is not one of its elements. In appendix 1, we provide a procedure 
to test the convex hull inclusion of a sensor, i.e., for any sensor, /, to determine if it lies in the convex hull of to + 1 
nodes arbitrarily chosen from the set, /C(/,r(), of its neighbors. Finding such a set O; is an important step in 
DILOC and we refer to it as triangulation and O; is referred to as a triangulation set. 

Let ci be the TO,-dimensional coordinates for a node, / G 6, with respect to a global coordinate system, written 
as the TO-dimensional row vector, 

C; = [C;,1,C;,2, • • • ,C;,to] . (4) 

The true (possibly unknown) location of sensor / is represented by . Because the distributed localization algorithm 
DILOC is iterative, c; {t) will represent the location vector, or state, for sensor I at iteration t. 

Barycentric coordinates. DILOC is expressed in terms of the barycentric coordinates, oj^, of a sensor, I e ^1, 
with respect to the nodes, k € 6;, see [2], [3]. The barycentric coordinates, aik, are unique and are given by 

aik = 1 , (5) 

'The convex hull, C ('J'), of a set of points in * is the minimal convex set containing ^. 
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with ^ 0, where '\' denotes the set difference, A^i}u@i\{k} is the generalized volume of the set {/} U6; \ {k}, 
i.e., the set 6; with sensor I added and node k removed. The barycentric coordinates can be computed from the 
inter-sensor distances dik using the Cayley-Menger determinants as shown in appendix II. From (5), and the facts 
that the generalized volumes are non-negative and 

Aeiu{l}\{k} = ^ € ^(6;), (6) 

fcee; 

it follows that, for each I G fl, k G Qi, 

aik e [0, 1], = 1- C^) 

B. Distributed iterative localization algorithm. 

Before presenting DILOC, we state expUcitly state our assumptions. 
(BO) Nondegeneracy. The generaUzed volume for k, A^^Q? 

(Bl) Anchor nodes. The anchors' locations are known, i.e., their state remains constant 

Cg(i) = c*, g e t > 0. (8) 
(B2) Convexity. All the sensors lie inside the convex hull of the anchors 

C(f2) C C{k). (9) 

From (B2), the next Lemma follows easily. 

Lemma 1 (Triangulation) For every I e fl, there exist r; > and 6; (r;) with |6( (r;)| = m+ 1 satisfying the 
properties in (3). 

Proof: Clearly, by (B2), k satisfies (3) and the diameter of the network, maxj^^ dik, {I £ il,k £ k), could be 
taken as r;. ■ 
Lemma 1 provides an existence proof, but in localization in wireless sensor networks, it is important to triangulate 
a sensor not with the network diameter but with a small r;. In fact. Section II-E discusses the probability of finding 
one such O; with r/ ^ maxjjt dik, G fi, fc e k). 

As a note, it is easily verified that every pair of sensors in C (O; (r;)) is within a distance Ri — 2r;. Think 
of m = 2, then |0( = 3, C (O; (r;)) is a triangle and if ri is the maximum distance of any interior point I to 
the vertices of the triangle, the maximum distance between the 4 points (triangle vertices and is 2ri. We complete 
stating the assumptions underlying DILOC. 

^Nondegeneracy simply states that the anchors do not He on a hyperplane. If this was the case, then the locaHzation problem reduces to 
a lower dimensional problem, i.e., R™"^ instead of R™. For instance, if all the anchors in a network lie on a plane in R^, the localization 
problem can be thought of as localization in R^. 
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(B3) inter-sensor distances. For I e fl, there exists at least an r; > and 9; (r;) C K{l,ri), satisfying (3), 
such that / has a communication link for every k G Qi (r;) and knows the mutual distances among all the nodes 
m{i}uei in). 

We can now present DILOC. There are two steps: a set-up phase and then DILOC proper. We discuss each 
separately. 

DILOC set-up: Triangulation. In the set-up phase, each sensor I triangulates itself, so that by the end of this 
phase we have paired every I e fl with its corresponding m+1 neighbors in 6;. Since triangulation should be with 
a small r;, the following is a practical protocol for the set-up phase. A sensor starts with a small conmiunication 
radius, r;, and chooses arbitrarily m+1 sensors within r; and tests if it hes in the convex hull of these sensors using 
the procedure described in Appendix I. The sensor attempts this with all collections of m -|- 1 sensors within ri. 
If all attempts fail, the sensor adaptively increases, in small increments, its communication radius, r;, and repeats 
the process. By (B2) and (9), success is eventually achieved and each sensor is triangulated by finding 6; with 
properties (3) and (B3). If the sensors have directionality a much simpler algorithm based on Lemma 2 below, see 
also discussion following the Lemma, triangulates the sensor with high probability of success in one shot. To assess 
the practical impUcations required by DILOC s set-up phase. Subsection II-E considers the realistic scenario where 
the sensors are deployed using a random Poisson distribution and computes in terms of deployment parameters the 
probability of finding at least one such 6; in a given radius, r;. 

DILOC iterations: state updating. Once the set-up phase is complete; at time t+1, each sensor / e f2, iteratively 
updates its state, i.e., its current location estimate, by a convex combination of the states at time t of the nodes 
in 6;. The anchors do not update their state, since they know their locations. The updating is explicitly given by 



where aik are the barycentric coordinates of I with respect to k G 6/. DILOC in (10) is distributed since (i) the 
update is implemented at each sensor independently; (ii) at sensor / G fi, the update of the state, ci{t +1), is 
obtained from the states of its m -|- 1 neighboring nodes in 8;; and (iii) there is no central location and only local 
information is available. 

DILOC: Matrix format. For compactness of notation, we write DILOC (10) in matrix form. Without loss of 
generality, we index the anchors in k as 1 , 2, . . . , m + 1 and the sensors in il as m + 2, m + 3, . . . , m + 1 + M = N. 
We now stack the (row vectors) states, C;, given in (4) for all the N nodes in the network in the N x m-dimensional 
coordinate matrix 




C;(i), I e K, 



(10) 



c = 



T T 
Ci , • • • , 



(11) 



DILOC equations in (10), now become in compact matrix form 



c(i-hi) = rc{t). 



(12) 
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Fig. 1. Deployment corresponding to tlie example in Section Il-C. 



The Structure of the N x N iteration matrix T becomes more apparent if we partition it as 



I, 



m+l 







T 



(13) 



B 



P 



The first m+l rows correspond to the update equations for the anchors in k. Since the states of the anchors are 
constant, see (B 1) and (8), the first m + l rows of T are zero except for a 1 at their {q,q),q e k = {l,...,r7i+l} 
location. In other words, these first m+l rows are the (m+l) x N block matrix [lm+i|0], i.e., the (m+l) x (m+l) 
identity matrix I„i+i concatenated with the (m + 1) x M zero matrix, 0. 

Each of the M remaining rows in T, indexed hy I E = {m + 2, m + 3, . . . , N}, have only m+l non- 
zero elements corresponding to the nodes in the triangulation set, Qi, of /, and these non-zero elements are the 
barycentric coordinates, aik, of sensor I with respect to the nodes in O;. The M x (m + 1) block B = {bij} is a 
zero matrix, except in those rows corresponding to the sensors in il that have a direct link to anchors. The M x M 
block P = {pij} is also a sparse matrix where the non-zero entries in row I correspond to the non-anchor nodes 
in 0;. The matrices T, P, and B have important properties that will be used to prove the convergence of the 
distributed iterative algorithm DILOC in Sections III and IV. 

Remark. Equation (12) writes DILOC in matrix format for compactness; it should not be confused with a 
centralized algorithm-it still is a distributed iterative algorithm. It is iterative, because each iteration through (12) 
simply updates the (matrix of) state(s) from C{t) to C(t + 1). It is distributed because each row equation updates 
the state of sensor I from the states of the m+l nodes in O/. In all, the iteration matrix, T, is highly sparse having 
exactly {m + 1) + A/(m + 1) non-zeros out of possible {m + 1 + M)^ elements. 



We consider an m = 2-dimensional Euclidean plane with m + 1 = 3 anchors and M = 4 sensors, see Fig. 1. 
The nodes are indexed such that the anchor set is k = {1,2,3}, |k| = m + 1 = 3, and the sensor set is f2 = 
{4, 5, 6, 7}, \n\ = M = 4. The set of all the nodes in the network is, thus, e = KUfi = {l,...,7}, |e| = iV = 7. 
The triangulation set, O;, I e il, identified by using the convex hull inclusion test are O4 = {1,5,7}, 65 = 



C. Example 
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{4,6,7}, Ge = {2,5,7}, 67 = {3,4,6}. It is clear that these triangulation sets satisfy properties in (3). It can 
be verified that sensor 5 does not have any anchor node in its triangulation set, 65, and every other sensor has 
only one anchor in its respective triangulation set. Since, no sensor is able to communicate to m + 1 = 3 anchors 
directly, no sensor can localize itself in a single step. 

At each sensor, I e fl, the barycentric coordinates, aik, k G Qi, are computed using the inter-sensor distances 
(among the nodes in the set {1} U 6;) in the Cayley-Menger determinant. It is noteworthy that the inter-sensor 
distances that need to be known at each sensor / to compute aik are only the inter-sensor distances among 
the m -|- 2 sensors in the set {1} U 6;. For instance, the distances in the Cayley-Menger determinant needed 
by sensor 5 to compute 054,055,057 are the inter-sensor distances among the nodes in the set, {5} U 65, i.e., 
c^54j £^565 dsT, d4Q, (^47, ^67. Duc to (3) and the relation 65 C /C(Z, r^), the nodes in {5} U65 he in a circle of radius, 
r5, centered at sensor 5 (shown in Fig. 1); on the other hand, no two sensors in the set {5} U 65 can be more than 
^5 = 2r5 apart. This justifies the choices of Ri and r; = Ri/2 in (2). 

Once the barycentric coordinates, aik, are computed, DILOC for the sensors in f2 is 

ci{t+l) = Y^aikCkit), /en = {4,5,6,7}. (14) 
keSi 

In particular, for sensor 5, we have the following expression, 

C5(i + 1) = a54C4(i) + a56C6(i) + a57C7(i)- 
DILOC for the anchors is given by 

Cg(t)=c*, qGK. 
We write DILOC for this example in the matrix format (12). 





-1) " 




1 
























C2(tH 


-1) 







1 



















C2(i) 


C3(n 


-1) 










1 
















caW 


C4(t H 


-1) 




041 











045 





047 




C4(t) 


C5(n 


-1) 













054 





056 


057 






C6(n 


-1) 







062 








065 





067 






C7(tH 


. 










073 


074 





076 










where the initial condition are C(0) = [0^,02,03,04,05,05,0°]^, with 0°, / £ Q, being randomly chosen row 
vectors of appropriate dimensions. 

Note here again that (15) is just a matrix representation of (14). DILOC is implemented in a distributed fashion 
as in (14). The matrix representation in (15) is included for compaction of notation and for the convergence analysis 
of the algorithm. 
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D. Complexity of DILOC 

Once the barycentric coordinates are computed, each sensor-^ performs the update in (10) that requires m + 1 
multipHcations and m additions. Assuming the computation complexity of the multiplication and the addition 
operations to be the same, the computation complexity of DILOC is 2m + 1 operations, i.e., 0(1) per sensor, per 
iteration. Since each sensor exchanges information with m + 1 nodes in its triangulation set, the communication 
complexity of DILOC is m + 1 communications, i.e., 0(1) per sensor, per iteration. Hence, both the computation 
and communication complexity are 0{M) for a network of M sensors. Note that, since the triangulation set-up 
phase, which identifies Qi (r;) and computes the barycentric coordinates, as explained in Subsection II-B, are to 
be carried out only once at the start of DILOC, they require a constant computation/communication complexity, so 
we do not account explicitly for it. 

E. Random Poisson Deployment 

A common model in wireless sensor networks is the Poisson deployment [27], [28]. We illustrate it on the plane, 
m — 2; the discussion can be extended to arbitrary dimensions. For a Poisson distribution with density, 7 > 0, 
the mean number of sensors in a sector Q with area Aq is 7v4q. The numbers of sensors in any two disjoint 
sectors, Qi and Q2, are independent random variables, and the locations of the sensors in a sector Q are uniformly 
distributed. We now characterize the probability of triangulating a sensor I in a region of radius, rj, centered at /. 
To this end, consider Fig. 2(a), which shows the triangulation region, a circle of radius, rj, centered at sensor /. Let 




(a) (b) 



Fig. 2. (a) Sensor / identifies its triangulation set, 0;, in the circle of radius, r;, centered at sensor I. The circle is divided into four disjoint 
sectors with equal areas, Qi , . . . , Q4. A suiScient condition for triangulation is that there exists at least one sensor in each of these four sectors, 
(b) Illustration of Lemma 2. 

Qi, Q2, Qs, Q4 be four disjoint sectors partitioning this circle with equal areas, i.e., Aq. — i — 1, . . . , 4. 

Lemma 2 A sufficient condition to triangulate a sensor Z e is to have at least one sensor in each of the four 
disjoint equal area sectors, Qi, i — 1, ... ,4, which partition the circle of radius of radius, r/, centered at I. 

^By sensor, we usually mean a non-anchor node. 



August 7, 2008 



DRAFT 



11 



Proof: In Fig. 2(b) consider the triangulation of sensor / located at the center of the circle; we choose arbitrarily 
four sensors Pi,P2,P3,P4 in each of the four sectors Qi, Q2, Q3, Q4. Denote the polygon with vertices Pi,P2,P3,P4 
by Pol{piP2P3P4). Consider the diagonal^ pi — p3 that partitions this polygon into two triangles Apip2P3 and 
ApiPaPi- Since / e Pol {P1P2P3P4.) and ApiP2P3 U Apip^p^ = Pol {P1P2P3P4) with ApiP2P3 n Apip^p^ = 0, 
then either I G Apip2P3 or I G Apip3P4. The triangle in which / Ues becomes the triangulating set, 6;, of I. 

This completes the proof. The generahzation to higher dimensions is straightforward; for instance, in R^, we 
have eight sectors and an arbitrary sensor I is triangulated with at least one sensor in each of these eight sectors 
(with equal volume) of a sphere of radius, ri, centered around I. ■ 

Let be the set of sensors in the sector Qi. It follows from the Poisson deployment that the probabihty of 
finding at least one sensor in a sector, Qi, of area Aq. is 

P(|QJ >0) = (l-exp-^^«^) . (16) 

Since the distribution of the number of sensors in disjoint sectors is independent, the probabihty of finding at least 
one sensor in each of the sets, Qi, . . . ,Q^, is the product 

P {\Qi\ > 0, Vi) = (1 - exp-^"'-' . (17) 

Thus, we have 

Pe, = P (6; exists satisfying (3)) > P (IQJ > 0, V i) . (18) 

This shows that, for a given deployment density, 7, we can choose, r;, appropriately, to guarantee the triangulation 
with arbitrarily high probabihty. Indeed, it follows from (18) that, for an arbitrary, < e < 1, to get the probability 
of triangulation to be greater than e, i.e., Pe, > e, the radius r; should be 

(19) 

In alternative, if we are limited by the communication radius, Ri, to guarantee Pe, > e, we will need the deployment 
density, 7, to be larger than 

7> ^^Infl-ei). (20) 

For example, if the sensors are deployed (following a Poisson distribution) with a density of 7 = 1 sensor/m^, 
then we can compute from equations (19)-(20) that 99% of the sensors will be able to triangulate (identify 6;) 
themselves if they can communicate to at least a radius of Ri = 5.52m. The remaining (1%) of the sensors may 
require to communicate to a larger radius. 

'*If Pol (P1P2P3P4) is concave, we choose the diagonal that lies inside Pol (pip2P3P4), i.e., pi — p3 in Fig. 2(b). If Pol(pip2P3P4) is 
convex, we can choose any of the two diagonals and the proof follows. 




August 7, 2008 



DRAFT 



12 



It also follows from the above discussion that if the sensors are equipped with a sense of directionahty (for 
example, if they all have ultrasound transducers) then each sensor has to find one neighbor in each of its 4 sectors, 
Qi,i,Qi,2, Qi,3, Qi,4 (in m = 2-d space). Once a neighbor is found, triangulation chooses 3 out of these 4, in order 
to identify 6;. The computational complexity in m = 2-d Euchdean space is 4 choose 3 = 4. Without directionality 
the process of finding 6; has the (expected) computation complexity of 'ynrf choose 3. 

III. Convergence of DILOC 

In this section, we prove the convergence of DILOC to the exact locations of the sensors, c*, I Gfl.To formally 
state the convergence result, we provide briefly some background and additional results, based on assumptions 
(B0)-(B3). 

The entries of the rows of the iteration matrix T, in (12), are either zero or the barycentric coordinates, aik, which 
are positive, and, by (7), add to 1. This matrix can then be interpreted as the transition matrix of a Markov chain. 
We then describe locahzation problem and DILOC in terms of a Markov chain. Let the assumptions (B0)-(B3) in 
Section II-B hold and the N nodes in the sensor network correspond to the states of a Markov chain where let the 
(y)-th element of the iteration matrix, T = {vij} defines the probabiUty that the ith state goes to the jth state. 
Because of the structure of T, this chain is a very special Markov chain. 

Absorbing Markov chain. Let an N x N matrix, T = {vij}, denote the transition probability matrix of a 
Markov chain with N states, Si,i=i,...,Ar. A state Sj is called absorbing if the probabihty of leaving that state is 
(i.e., Vij = 0,i ^ j, in other words vu = 1). A Markov chain is said to be absorbing if it has at least one 
absorbing state, and if from every state it is possible to go with non-zero probability to an absorbing state (not 
necessarily in one step). In an absorbing Markov chain, a state which is not absorbing is called transient. For 
additional background, see, for example, [29]. 

Lemma 3 The underlying Markov chain with the transition probability matrix given by the iteration matrix, T, is 
absorbing. 

Proof: We prove by contradiction. Since Wjj = 1, i e k, the anchors are the absorbing states of the Markov 
chain. Since vu = 0,i € O, the (non-anchor) sensors are the transient states. Partition the transient states into 
two clusters CI and C2, such that each transient state in CI can go with non-zero probability to at least one of 
the absorbing states and, with probabihty 1, the transient states in C2 cannot reach an absorbing states. It foUows 
that with probability 1 the transient states in C2 cannot reach the transient states in CI (in one or multiple steps); 
otherwise, they reach an absorbing state with a non-zero probabihty. Let's consider the he on the boundary of the 
convex huU, C(C2), i.e., the vertices of C(C2). Because they are on boundary, they cannot lie in the interior of the 
convex huU any subset of sensors in C(C2), and, thus, cannot triangulate themselves, which contradicts Lemma 1 
and assumption (B3). In order to triangulate the boundary sensors in C(C2), the boundary sensors in C2 must be 
able to reach the transient states and/or the absorbing states, that is to say that the boimdary sensors in C(C2) have 
to reach the sensors in CI to be able to triangulate themselves. Hence, the Markov chain is absorbing. ■ 
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Consider the partitioning of the iteration matrix, T, in (13). With the Markov chain interpretation, the M x (m+ 1) 
block B = {bij} is a transition probability matrix for the transient states to reach the absorbing states in one-step, 
and the block M x M P = {pij} is a transition probabiUty matrix for the transient states. With (13), T*+^ can be 
written as 



-i+l 



l-m+1 







fe=0 



(21) 



and, as t goes to infinity, we have 



lim T'+^ 



(22) 



Im+1 

(Im-P)"'B 

by Lemmas 5 and 6, in appendix HI. Lemmas 5 and 6 use the fact that if P is the matrix associated to the transient 
states of an absorbing Markov chain, then p{P) < 1, where p{-) is the spectral radius of a matrix. With (22), 
DILOC (10) converges to 

Im+1 

(Im-P)-'B 

In (23), the coordinates of the M sensors in il (last M rows of C(t + 1)) are written as a function of the m + 1 
anchors in k whose coordinates are exactly known. The limiting values of the states of the M sensors in are 
written in terms of the coordinates of the m+ 1 anchors in k weighted by (Im — P)~^B. To show that the limiting 
values are indeed the exact solution, we give the following Lemma. 



lim C{t + 1) 

t^OO 



C(0). 



(23) 



Lemma 4 Let c* be the exact coordinates of a node, I G 6. Let the M x (m+l) matrix, D = {dij}, I G G k, 
be the matrix of the barycentric coordinates of the M sensors (in Cl) in terms of the m + 1 anchors in k, relating 
the coordinates of the sensors to the coordinates of the anchors by 



I e n. 



(24) 



Then, we have 



D = (Im-P)"'B. 



(25) 



Proof: Clearly {1m — P) is invertible, since, by (95) in Appendix 111, p{P) < 1; this follows from the 
fact that the eigenvalues of the matrix Im — P are 1 — Xj, where is the jth eigenvalue of the matrix P and 
I Aj I < 1, Vj = 1, . . . , M. It suffices to show that. 



D = B + PD, 



(26) 



since (25) follows from (26). In (26), D and B are both M x (m + 1) matrices, whereas P is an M x M matrix 
whose non-zero elements are the barycentric coordinates for the sensors in Hence, for the Ij-th element in (26), 
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we need to show that 



dij = bij + y^^pikdkj- 



(27) 



For an arbitrary sensor, I e Q, its triangulation set, 6;, may contain nodes from both k and Q. We denote Ke, 
as the elements of 6; that are anchors, and as the elements of 6; that are non-anchor sensors. The exact 
coordinates, c*, of the sensor, I, can be expressed as a convex combination of the coordinates of its neighbors in 
its triangulation set, k G 9;, using the barycentric coordinates, aik, i.e., 

ke&i 

= YbijC* + J2Pikcl (28) 
j&K fceo 

since the scalars, aij, are given by 



bij, ifje/te,, 
o-ij = <t Pij, if i e fie,, 

0, ifj^e^. 

Equation (28) becomes, after writing each k Gfl in terms of the m + 1 anchors in k. 



(29) 



(30) 



This is a representation of the coordinates of sensor, I, in terms of the coordinates of the anchors, j G k. Since for 
each j e K, the value inside the parentheses is non-negative with their sum over j G k being 1 and the fact that 
the barycentric representation is unique, we must have 



dij = bij + y^^pikdkj, 
keo 

which, comparing to (24), completes the proof. 

We now recapitulate these results in the following theorem. 



(31) 



Theorem 1 (DILOC convergence) DILOC (10) converges to the exact coordinates, c*, of the M sensors (with 
unknown locations) in CI, i.e., 

]im ci{t + 1) = c* , VZefi. (32) 

t — ^oo 
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Proof: The proof is a consequence of Lemmas 3 and 4. ■ 
Convergence rate. The convergence rate of the locaUzation algorithm depends on the spectral radius of the 
matrix P, which by (95) in Appendix III is strictly less than one. This is a consequence of the fact that P is a 
sub- stochastic matrix. The convergence is slow if the spectral radius, p(P), is close to 1. This can happen if the 
matrix B is close to a zero matrix, 0. This is the case if and only if the sensors cluster in a region of very small area 
inside the convex hull of the anchors, and the anchors themselves are very far apart. In fact, it can be seen that in 
this case the barycentric coordinates for the sensors with kq, ^ (see Lemma 4 for this notation) corresponding 
to the elements in kq, are close to zero. Since in practical wireless sensor applications the sensors are assumed to 
be deployed in a geometric or a Poisson fashion (see details in Section II-E), the probabiUty of this to happen is 
arbitrarily close to 0. 



IV. DILOC WITH Relaxation 

In this Section, we modify DILOC to speed its convergence rate and to obtain a form that is more suitable to 
study distributed localization in random environments. We observe that in DILOC (10), at time t+l,ihe expression 
for ci{t + 1), I G Cl, does not involve its own coordinates, ci{t), at time t. We introduce a relaxation parameter, 
a e (0, 1], in the iterations, such that, the expression of ci{t + 1) is a convex combination of ci{t) and (10). We 
refer to this version as the DILOC with relaxation, DILOC-REL. It is given by 



ciit+1) = 



{1 - a)ci{t) + aci{t) ^ ci{t), I e K, 



(1 - a)ci{t) +aJ2kee, aikCk{t), I G Q.. 
DILOC is the special case of DILOC-REL with a = 1. The matrix representation of DILOC-REL is 



(33) 



C(i + 1) =HC(t), 



(34) 



where H = (1 — a) Ijv + aX and Ijv is the N x N identity matrix. It is straightforward to show that the iteration 
matrix, H, corresponds to a transition probability matrix of an absorbing Markov chain, where the anchors are the 
absorbing states and the sensors are the transient states. Let J = (1 — a) Im + aP; partitioning H as 



We note the following 



and, as t — > oo. 



H 



H 



Im+1 
QfB J 

Im+1 



(35) 



k=0 



(36) 



lim H*+i 



Im+1 

(lM-J)"^aB 



(37) 
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from Lemmas 5 and 6. Lemmas 5 and 6 apply to H, since H is non-negative and p(J) < 1. To show p(J) < 1, 
we recall that p(P) < 1 and the eigenvalues of J are (1 — a) + aXj, where Xj are the eigenvalues of P. Therefore, 
we have 

p(J) = max|(l - a) +aAj| < 1. (38) 
j 

The following Theorem establishes convergence of DILOC-REL. 

Theorem 2 DILOC-REL (33) converges to the exact coordinates, c*, of the M sensors (with unknown locations) 
in CI, i.e., 

]im ci{t + 1) = c* , VZefi. (39) 

Proof: It suffices to show that 

(lM-J)"'aB = (Im-P)"'B. (40) 

To this end, we note that 

(lM-J)"'aB = (lM-((l-a)lM + aP))-'aB, (41) 

which reduces to (40) after basic algebraic manipulations. The convergence of DILOC-REL thus follows from 
Lemma 4. ■ 
As mentioned, the advantage of DILOC-REL is twofold: since p{3) is a function of a, we may optimize the 
convergence rate over a; and DILOC-REL forms the basis for the distributed localization algorithm in random 
environments (DLRE) that we discuss in Sections V and VI. 

V. Distributed Localization in Random Environments: Assumptions and Algorithm 

This and the next Section study distributed iterative localization in more realistic practical scenarios, when the 
inter-sensor distances are known up to errors, the communication links between sensors may fail and, when alive, 
the communication among sensors is corrupted by noise. We write the update equations for DILOC-REL, (34), in 
terms of the columns, c^(t), 1 < j < m, of the coordinate matrix, C(t). Column j corresponds to the vector of 
the j-th estimate coordinates of all the N sensor locations^. The updates are 

c^{t + l) = [{l-a)I + ar] c> {t), l<j< m. (42) 

We partition c^{t) as 

c'{t)= , (43) 

^In the sequel, we omit the subscripts from the identity matrix, I, and its dimensions will be clear from the context. 
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where, u-' e corresponds to the j'-th coordinates of the anchors, which are know (hence, we omit the 

time index, as they are not updated), and x-'(f) e M^^^ corresponds to the estimates of the j-th coordinates of 
the non-anchor sensors, hence not known. Since, update is performed only on the x^{t), (42) is equivalent to the 
following recursion: 

x^{t+l) = [{l-a)I + aP] (i) + aBu^ . (44) 

Thus, to implement the sequence of iterations in (44) perfectly, the l-th sensor at iteration t needs the corresponding 
rows of the matrices P and B, and, in addition, the current estimates, c^(t),n G G; (j-th component of the n-th 
sensor coordinates), of its neighbors' positions. In practice, there are several limitations: (i) The computation of 
the matrices P and B requires inter-sensor distance computations, which are not perfect in a random environment; 
(ii) the communication channels, or links, between neighboring channels may fail at random times; and (iii) because 
of imperfect communication, each sensor receives only noisy versions of its neighbors current state. Hence, in 
a random environment, we need to modify the iteration sequence in (44) to account for the partial imperfect 
information received by a sensor at each iteration. We start by stating formally our modeling assumptions. 
(CI) Randomness in system matrices. At each iteration, each sensor needs the corresponding row of the system 
matrices B and P, which in turn, depend on the inter-sensor distance measurements, which can be, possibly, 
random. Since a single measurement of the inter-sensor distances may lead to a large random noise, we assume 
the sensors estimate the required distances at each iteration of the algorithm (note that this leads to an implicit 
averaging of the unbiased noisy effects, as will be demonstrated later.) In other words, at each iteration, the l-th 
sensor can only get estimates, B;(t) and Pi{t), of the corresponding rows of the B and P matrices, respectively. 
In the generic imperfect commimication case, we have 



where 



{SbW| 
L J t>o 



B(t) =B + SB + SB(i) 
is an independent sequence of random matrices with. 



E 



SB(t) 



= 0, yt, supE 

t>0 



SbW 



= fcs < 00. 



Here, Sb is the mean measurement error. Similarly, for P, we have 

P(t) = P + Sp + Sp(t), 

where I Sp (t) > is an independent sequence of random matrices with, 
I. J t>o 



E 



Sp(t) = 0, \/t, supE 
t>o 



Sp(i) 



fcp < oo. 



(45) 



(46) 



(47) 



(48) 



Likewise, Sp is the mean measurement error. Note that this way of writing B{t), P{t) does not require the noise 
model to be additive. It only says that any random object may be written as the sum of a deterministic mean part 
and the corresponding zero mean random part. The moment assumptions in eqns. (46,48) are very weak and, in 
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particular, are satisfied if the sequences i B(t) > and {Fit) > are i.i.d. with finite variance. 

L J t>o L J t>o 

(C2) Random Link Failure: We assume that the inter-sensor communication links fail randomly. This happens, for 
example, in wireless sensor network applications, where occasionally data packets are dropped. To this end, if the 
sensors I and n share a communication Unk (or, n G 6;), we assume that the hnk fails with some probabihty 1 — qin 
at each iteration, where < qin < 1. We associate with each such potential network Unk, a binary random 
variable, ein{t), where e;„(t) = 1 indicates that the corresponding network link is active at time t, whereas e;„(f) = 
indicates a hnk failure. Note that E [e;„] = g;„. 

be a family of independent zero mean random variables such 



(C3) Additive Channel Noise: Let 

that 



sup E vlit) 



ky < oo. 



(49) 



We assume that, at the t-th iteration, if the network link {I, n) is active, sensor I receives only a corrupt version, 
yjn{t), of sensor n's state, c^(i), given by 



(50) 



This models the channel noise. The moment assumption in eqn. (49) is very weak and holds, in particular, if the 
channel noise is i.i.d. 

(C4) Independence: We assiune that the sequences, \ Ss{t),Sp{t) [ , {e(„(t)}; „ ^, and < vj^{t) \ are 
mutually independent. These assumptions do not put restrictions on the distributional form of the random errors, 
only that they obey some weak moment conditions. 
Clearly, under the random envirormient model (as detailed in Assumptions (C1)-(C4), the algorithm in (44) is not 
appropriate to update the sensors states. We now consider the following state update recursion for the random 
envirormient case. 

Distributed Localization in Random Environment Algorithm (DLRE): 

ein{t)Bin{t) 



xi{t+l) = {l-a{t))xj{t) + a{t) 



E 



E 



ein{t)hn{t) 



qin 



qin 



, I G^l, l< j <m 



(51) 



In contrast with DILOC-REL, in (51), the gain a{t) is now time varying. It will become clear why when we 
study convergence of this algorithm. To write DLRE in a compact form, we introduce notation. Define the random 
matrices, B{t) e M^x(™+i) and P(i) e M^xC^+i), as the matrices with In entries given by 



B;„(t) = B;„(t) 



einjt) 

qin 



1 , Pin{t) = Pln{t) 



einit) 



qir. 



(52) 
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Clearly, by (C2),(C4), the matrices B(i) e RMx{m+i) p(^) g ^Mx{m+i) ^ero mean. Note that E [e(„] = 
qin- Also, by the bounded moment assumptions in (CI), we have 



supE 
t>o 



= fcs < oo, supE 
t>o 



Fit) 



= fcp < oo. 



Hence, the iterations in (51) can be written in vector form as 

x^ (t + 1) = (1 - a{t)) x^ (i) + a{t) [(P(i) + P(t)) x^ (t) + (B(i) + B(i)) u^' + 7f{t) 
where, the Zth element of the vector, ri^{t), is given by 

By (C1)-(C4), the sequence, {??-' (i)}t>o. is zero mean, independent, with 

|2" 



(53) 



(54) 



(55) 



supE \\rf{t)\ 
t L 



— < OO. 



(56) 



(57) 



From (CI), the iteration sequence in (54) can be written as 

jc\t + 1) = x^ (t) - a{t) [(I - P - Sp) x-''(t) - (B + Sb) u^' - (Sp(i) + P(i)) x^ (t) 
-(SB(i)+B(t))u^"-77^(f)" . 

We now make two additional design assumptions. 
(Dl) Persistence Condition: The weight sequence satisfies 

a{t)>0, ^a{t)=oo, ^a^it) < OO. (58) 

t>0 i>0 

This condition, commonly assumed in the adaptive control and adaptive signal processing literature, assumes that 
the weights decay to zero, but not too fast. 
(D2) Low Error Bias: We assume that 

p(P + Sp)<l. (59) 

Clearly, we have p(P) < 1. Thus, if we assume that the non-zero bias, Sp, in the system matrix (resulting 
from incorrect distant computation) is small, (59) is justified. We note that this condition ensures that the matrix 
(I — P — Sp) is invertible. 

In the following sections, we prove that the DLRE algorithm, under the assumptions (C1)-(C4), (D1)-(D2), leads 
to a.s. convergence of the state vector sequence, {x^ (t)}^^g, to a deterministic vector for each j, which may be 
different from the exact sensor locations, because of the random errors in the iterations. We characterize this error 
and show that it depends on the non-zero biases, Sb, and Sp in the system matrix computations, and vanishes as 
IISbII ^ and ||Sp|| 0. 
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VI. DLRE: A.S. CONVERGENCE 
We show the ahnost sure convergence of DLRE under the random environment presented in Section V. 

Theorem 3 (DLRE a.s. convergence) Let {-xP {t)}t>o-i 1 < i < be the state sequence generated by the iterations, 
given by (57), under the assumptions (C1)-(C4), (D1)-(D2). Then, 

P [ lim x^ (t) = (I - P - Sp)"^ (B + Sb) Vjl = 1. (60) 

The convergence analysis of the DLRE algorithm is based on the sample path properties of controlled Markov 
processes, which has also been used recently to prove convergence properties of distributed iterative stochastic 
algorithms in sensor networks, e.g., [30], [31]. The proof relies on the following result from [32], which we state 
here as a theorem. 

Theorem 4 Consider the following recursive procedure: 

x(i + 1) = x(t) + a{t) [R(x(i)) + T{t + 1, x(i), w)] , (61) 

where, x, R, F are vectors in K^^^. There is an underlying conmion probability space (H,^, P), and oj is the 
canonical element of the probability space, H. Assume that the following conditions are satisfied^. 

1) : The vector function R(x) is Borel measurable and V{t, x, oj) is B'^ T measurable for every t. 

2) : There exists a filtration {Tt}t>o of such that the family of random vectors r(f,x, w) is measurable, 
zero-mean and independent of Tt-i- 

3) : There exists a function V{-x) e with bounded second order partial derivatives satisfying: 

y(xo) = 0, t/(x) > 0. x^xo. (62) 
sup (R(x),14(x)) < 0, Ve>0. (63) 

||x-xo||>e 

4) : There exist constants fci, ^2 > 0, such that, 

||R(x)||2 + E[||r(i,x,a;)||2] <fci(l + y(x))-fc2(R(x),y,,(x)) (64) 

5) : The weight sequence {a(t)}t>o satisfies the persistence condition (Dl) given by (58). 
Then the Markov process, {x(t)}t>o, converges a.s. to xq. 

Proof: The proof follows from Theorem 4.4.4 in [32] and is omitted due to space constraints. ■ 
We now return to the proof of Theorem 3. 

*In the sequel, B'^ denotes the Borel sigma algebra in R**^^!. The space of twice continuously differentiable functions is denoted by C2, 
while 14: (x) denotes the gradient ^ g^^^ . 
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Proof: [Proof of Theorem 3] We will show that, under the assumptions, the algorithm in (57) falls under the 
purview of Theorem 4. To this end, consider the filtration, {^t}t>o. where 



Define the vector d* as 



= a (x^(0),Sp(s),P(s),SB(s),B(s),7y^(s) : < s < . 



d* = (I - P - Sp)-' (B + Sb) u^'. 



Equation (57) can be written as 

jc\t+l)=^\t) - [(I-P-Sp)(x^(t)-d*)- (Sp(t) + P(f))x^(<) 

- (SBW + B(t))u^ -r?^"(t) . 

In the notation of Theorem 4, (67) is given by 

x^(t + l) =x^(t)+a(t) [R(x^(t))+r(t+l,x^(i),w)] , 

where 

R(x) = - (I - P - Sp) (xJ'(i) - d*) , 

and 

r(i + l,x,w)= [(Sp(i)+P(t))x^(t)+ (sBW + B(t))u^ +r?^(t) . 



(65) 



(66) 



(67) 



(68) 



(69) 



(70) 



This definition satisfies assumptions 1) and 2) of Theorem 4. 

We now show the existence of a stochastic potential function V{-) satisfying the remaining assumptions of 
Theorem 4. To this end, define 



y(x) = ||x-d* 



Clearly, V{x) e C2 with boimded second order partial derivatives, with 



(71) 



V{d*) = 0, F(x) >0, x^d*. 



(72) 
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sup (R(x),yx(x)) 

l|x-dJ|>e 



< 



< 



sup -2(x-d*)^(I-P-Sp)(x-d*), 

|x-d.||>e 



2 (x - d*)"^ (P + Sp) (x - d*) - 2||x - d* 



(x-d*)^(P + Sp)(x-d*) -2||x-d 



]*|l2 



sup 

|x-d.||>e 

sup 

|x-d.||>e 

sup [2 ||x - d*|| p(P + Sp) ||x - d*|| - 2||x - d*f] , 

lx-dJ|>e 



= sup -2(l-p(P + Sp))||x-d*f , 

l|x-d.||>e 

< -2e2(l-p(P + Sp)), 

< 0, 



(73) 



where, the last step follows from (D2). Thus, assumption 3) in Theorem 4 is satisfied. 
To verify 4) note that 

||R(x)f = (x-d*)^(I-P-Sp)^(I-P-Sp)(x-d*), 

< |(I-P-Sp)^(I-P-Sp) 

= A;i||x-d*||', 

= fciy(x), 



(74) 



where fci > is a constant. 

Finally, by assumptions (C1)-(C4), we have 



E||r(t,x,a;)f = E [(^Sp(t - 1) + P(t - 1)] x + {Ssit - 1) + B{t - l)j + rf {t - 1)\ 
[(Sp(t - 1) + P{t - 1)) X + (SbC* - 1) + B(t - 1)) u^' +v'{t- 1)] , 
= x'^E [Sp(t - l)Sp(t - 1) + P^{t - l)P(t - 1)] X + u-''^E [sS(t - l)SB(t - 1) 

+ B^(t- l)B(t- 1)] u^' +E ||r/^(t- +2x^E [Sp(t- l)SB(t- 1)] u-'', 



< E 

+E 



Sp(t-l) + P(t-l) 



iixr+E 



isB(t-i)ir+iiB(t-i)ii' 



Wit -I) 



- 2E 



Sp(f- 1) 



1/2 



SB(t- 1) 



1/2 



< (fcp + fcp) ||xf + (kB + kB^ ||u-'||' +fe^ + fcpA;B ||x|| ||u-' 



(75) 



The cross terms dropped in the second step of eqn. (75) have zero mean by the independence assumption (C4). 
For example, consider the term E Sp{t— l)P(t — 1) .It follows from eqns. (47,52) that the In-th entry of the 
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matrix Sp{t — l)P(i — 1) is given by 

[s?(t-i)P(t-i)] = ^[s?(^-i)l, [p(t-i)l 

r 



rn \ t/rn 



+ Epp(^-i)],[sp(t-i)], 



6rn (^) 



6rn (^) 



From the independence and zero-mean assumptions, we have the following, Vr: 

E[[Sp(t-l)]jPu(^-l)] = [pue[[sp(*-i)]Je[ 

= 

E[[Sp(t-l)]jSpu(^-l)] = [SpUE[[Sp(t-l)]jE[(^-l)] 

= 

E [[S., - [S., - = E [[S., - [S., - 1)] E [(^ - .)] 



= 



where we have repeatedly used the fact that 



E 



From eqns. (76-77) it is then clear that 



E S?(i- l)P(t- 1) =0 



In a similar way, it can be shown that the other dropped crossed terms in eqn. (75) are zero-mean. 
We note that there exist constants, kz,kA,k^,kQ > 0, such that 

||x||2 < fcsllx - d*||^ + ki, ||x|| < fcsllx - d*||^ + ka. 



Hence, from (75) and (79), we have 



E||r(t,x,a;)||^ < fc7||x-d*f -Ffcg, 



< fc9(i + y(x)), 



(76) 



(77) 



(78) 



(79) 



(80) 



where, k-j,k^ > and fcg = max (fcy, /cg). Combining eqns. (74,80) we note that assumption 4) in Theorem 4 is 
satisfied, as 

(R(x),K(x)) <0, Vx. (81) 
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Hence, all the conditions of Theorem 4 are satisfied and we conclude that 



lim x^{t) = (I - P - Sp)"^ (B + Sb) 



1. 



Since, (82) holds for all j, and j takes only a finite number of values (1 < j < m), we have 



lim x^{t) = (I - P - Sp)"' (B + Sb) Vj 



= 1. 



(82) 



(83) 



We now interpret Theorem 3. Referring to the partitioning of the C{t) matrix in (43), we have 



u 



(84) 



where each row of X(f) corresponds to an estimated sensor location at time t. Theorem 3 then states that, starting 
with any initial guess, X(0) G M*^^"*, of the unknown sensor locations, the state sequence, {'K{t)}t>o, generated 
by the DLRE algorithm converges a.s., i.e.. 



lim X{t) = (I - P - Sp)"^ (B + Sb) U 

t^OO 



1. 



However, it follows from Lenmia 4, that the exact locations of the unknown sensors are given by 

X* = (I-P)"^BU. 



(85) 



(86) 



Thus, the steady state estimate given by the DLRE algorithm is not exact, and, to characterize its performance, we 
introduce the following notion of locahzation error, ei, as 



ei = 



(I - P - Sp)-' (B + Sb) U - (I - P)-' BU 



(87) 



We note that ei is only a function of Sp, Sb, the non-zero biases in the system matrix computations, resulting from 
noisy inter-sensor distance measurements (see, Section V.) We note that the DLRE algorithm is robust to random 
Unk failures and additive channel noise in inter-sensor communication. In fact, it is also robust to the zero-mean 
random errors in the system matrix computations, and only affected by the fixed non-zero biases. Note that e; = 
for Sp = Sb = 0. Clearly, if we assume sufficient accuracy in the inter-sensor distance computation process, so 
that the biases, Sp,Sb, are small, the steady state error Bss wiU also be negligible, even in a random sensing 
envirormient. These are illustrated by numerical studies provided in Section VII. 

VII. Numerical Studies 

We divide the numerical study of DILOC into the following parts. First, we present DILOC in the deterministic 
case, i.e., we have no communication noise, no link failures, and the required inter-sensor distances are known 
precisely. Second, we present DILOC when there is conmiunication noise and link failures. Third, we consider 
noise on the distance measurements that results in random system matrices, P{t) and B(t), as given in (45)-(47); 
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(a) 



(b) 



(c) 



(d) 



Fig. 3. Deterministic environments: (a) and (b) DILOC algorithm implemented on the example in Section II-C. (c) An = 500 node network 
and the respective triangulation sets, (d) DILOC implemented on the network in (c), where the iterations are shown for two arbitrarily chosen 
sensors. 





500 1000 1500 2000 
DILOC iteration, * 



(a) 



(b) 



Fig. 4. Effect of communication noise and link failures: (a) An A'^ = 50 node network and the respective triangulation sets, (d) DLRE (with a 
decreasing weight sequence, a = implemented on the network in (a), where the iterations are shown for two arbitrarily chosen sensors. 



we study both the biased (Sp ^ 0, Sb ^ 0) and unbiased (Sp = Sb = 0) cases in the following. In the end, we 
present studies of DILOC in the presence of all random scenarios. 

DILOC Algorithm in Deterministic Environments: We consider the example presented in section II-C. We 
have N — 7 nodes in m = 2-dimensional space, where m + 1 = 3 are the anchors and M — i are the sensors. 
DILOC, as given in (14), is implemented, and the results are shown in Fig. 3(a)-3(b). Fig. 3(a) shows the estimated 
coordinates of sensor 6 and Fig. 3(b) shows the trajectories of the estimated coordinates for all the sensors with 
random initial condition. Fig. 3(c) and Fig. 3(d) show DILOC for a network of = 500 nodes. 

DILOC Algorithm with Communication Noise and Link Failures: We consider the same example of section II- 
C, but, include noise in the communication and link failures. All the communication links are active 90% of the 
time, i.e., qin — 0.9, V n s.t. n G Qi, as discussed in (C2), and include an additive communication noise that is 
Gaussian i.i.d. with zero-mean and variance 1/M (roughly speaking, this is equivalent to having a unity variance 
for the entire network). In this scenario, we employ DILOC with a decreasing weight sequence, a{t) = and 
the results are presented in Fig. 4(a)and Fig. 4(b). 

DILOC with Noisy Distance measurements: We now consider noise on the distance measurements. We assume 
that we have a reasonable estimate of the required distances such that it translates into a small perturbation of the 
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(a) 



(b) 



Fig. 5. Effect of noisy distance measurements: (a) An A'^ = 50 node network and the respective triangulation sets, (d) DLRE (witli a decreasing 
weight sequence, a = j^^) implemented on the networlc in (a), where the iterations are shown for two arbitrarily chosen sensors. 





(a) 



(b) 



Fig. 6. Random environments (Noisy distances, communication noise, link failures): (a) An Af = 50 node network and the respective 
triangulation sets, (d) DLRE (with a decreasing weight sequence, a = ^trgg) implemented on the network in (a), where the iterations are 
shown for two arbitraiily chosen sensors. 



system matrices, B(t) and P(i). The matrices Sp,Sb are zero-mean Gaussian i.i.d. perturbations with variance 0.1 
(small signal perturbation, note that the non-zero elements of both the matrices, B and P lie in the range of 
and 1). For a network of = 50 nodes in Fig. 5(a), we implement DLRE, with a decreasing weight sequence, a = 
^0^, in Fig. 5(b). Finally, Fig. 6(a) shows a network of = 50 nodes, where DLRE, with a decreasing weight 
sequence, a = ^o^j with all of the above random scenarios is implemented in Fig. 6(b). 



VIII. Conclusions 

The paper studies a distributed iterative sensor localization algorithm in m— dimensional Euclidean space, M™ (m > 
1), that finds the location coordinates of the sensors in a sensor network with only local communication. The 
algorithm uses the minimal number, m + 1, of anchors (sensors with known location) to localize an arbitrary 
number, M, of sensors that lie in the convex hull of these m+1 anchors. In the deterministic case, i.e., when no noise 
affects the inter-sensor communication, the inter-sensor distances are known with no errors, and the communication 
links do not fail, we show that our distributed algorithms, DILOC and DILOC-REL, lead to convergence to the 
exact sensor locations. For the random environment scenario, where inter-sensor communication links may fail 
randomly, transmitted data is distorted by noise, and inter-sensor distance information is imprecise, we show that 
our modified algorithm, DLRE, leads to almost sure convergence of the iterative location estimates, and in this 
case we explicitly characterize the resulting error between the exact sensor locations and the converged estimates. 
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(a) 



(b) 



Fig. 7. Convex Hull Inclusion Test (m=3): The sensor I is shown by a 'o', whereas, the anchors in k are shown by 'V. (a) I S C{k) 



Numerical simulations illustrate the behavior of the algorithms under different field conditions. 

Appendix I 
Convex Hull Inclusion test 

We now give an algorithm that tests if a given sensor, I £ M'", lies in the convex hull of m + 1 nodes in a set, k, 
using only the mutual distance information among these m + 2 nodes (k U {I}). Let k denote the set of to + 1 
nodes and let C{k) denote the convex hull formed by the nodes in k. Clearly, if I E C{k), then the convex hull 
formed by the nodes in k is the same as the convex hull formed by the nodes in k U {/}, i.e.. 



C{k)^C{kU{1}), if I eC{K). 



(88) 



With the above equation, we can see that, if / e C(k), then the generalized volumes of the two convex sets, C{k) 
and C{kU{1}), should be equal. Let denote the generalized volume of C{k) and let denote the generalized 

volume of C{k U {I}), we have 



Hence, the test becomes 



= ^KU{i}, 



if I e C{k). 



I e C(k), if X! ^'^u{;}\{fc} 
I i C(k), if ^A^^{i}\{k} > A 

kEK 



(89) 

(90) 
(91) 



This is also shown in Figure 7. The above inclusion test is based entirely on the generalized volumes, which can 
be calculated using only the distance information in the Cayley-Menger determinants. 
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Appendix II 
Cayley-Menger Determinant 

Let K be a set of m + 1 points (sensors) in M™, and dy be the inter-sensor distance between I and j. The 
generalized volume, A^, of the convex hull of the points in k can be computed by the Cayley-Menger determinant, 
see, e.g., [33]. The Cayley-Menger determinant is the determinant ofanm-|-2xm-|-2 (synmietric) matrix that 
relates to the generahzed volume, A^, of the convex huU, C(k;), of the m -|- 1 points in through an integer 
sequence, Sm+i ■ The Cayley-Menger determinant is given by 






1 


1 


1 






1 


1 













r72 

"l,m+l 


1 
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"23 


"24 




r72 

"2,m+l 
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dh 
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dh 







^72 

"m-l,m+l 














"m,m+l 


1 


"m+1,1 


"m+1,1 




"TO+l,m-l 


^72 

"m+l,m 







= (-l)'"+^2™(m 




m = 


{0,1,2,.. 





(92) 



where 

(93) 

and its first few terms are -1, 2, -16, 288, -9216, 460800, .... 

Appendix 111 
Important Results 

Lemma 5 If the matrix, P, corresponds to the transition probability matrix associated to the transient states of an 
absorbing Markov chain, then 



Proof: For such a matrix, P, we have 



lim P*+i = 0. 



P(P) < 1, 



(94) 



(95) 



from Lemma 8.3.20 and Theorem 8.3.21 in [34], where p{-) denotes the spectral norm of a matrix and (94) follows 
from (95). ■ 

Lemma 6 If the matrix, P, corresponds to the transition probability matrix associated to the transient states of an 
absorbing Markov chain, then 

t+i 

limVp'^ = (I-P)"\ (96) 

/s=0 



Proof: 
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The proof follows from Lemma 5 and Lemma 6.2.1 in [34]. ■ 
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